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We construct phenomenologically viable supersymmetric models where CP is an approx- 
imate symmetry. The full high energy theory has exact CP and horizontal symmetries 
that are spontaneously broken with a naturally induced hierarchy of scales, Acp <C A#. 
Consequently, the effective low energy theory, that is the supersymmetric Standard Model, 
has CP broken explicitly but by a small parameter. The Ek parameter is accounted for 
by supersymmetric contributions. The predictions for other CP violating observables are 
very different from the Standard Model. In particular, CP violating effects in neutral B 
decays into final CP eigenstates such as B — > ijjKs and in K — > -kvv decays are very small. 



1. Introduction 



Within the Standard Model, the following features regarding CP violation hold: 

(i) CP is broken explicitly. 

(ii) All CP violation arises from a single phase (that is the Kobayashi-Maskawa phase 
5km)- 

(hi) The measured value of Sk requires that 5km is of order one. (In other words, CP is 

not an approximate symmetry of the Standard Model.) 
(iv) The values of all other CP violating observables can be predicted. In particular, the 

asymmetry a^K S t 

a, Ks S ,n(Am B *) - - - — - - ^ , (1.1) 

(and similarly various other CP asymmetries in B decays), and the ratio a nvj >, 

T(K L -> ir°uv) 



;i.2) 



are expected to be of order one. 

The commonly repeated statement that CP violation is one of the least tested aspects 
of the Standard Model is well demonstrated by the fact that none of the above features 
necessarily holds in the presence of New Physics. Such a dramatic difference from the 
Standard Model is possible, for example, in the supersymmetric framework. (For a recent 
review of CP violation in supersymmetry, see [1J.) Indeed, in this work, we construct 
phenomenologically viable supersymmetric models, with the following features: 

(i) CP is an exact symmetry of the full theory but is spontaneously broken at some high 
energy scale by a VEV of a gauge singlet scalar field. 

(ii) In the low energy effective theory, there are many independent CP violating phases, 
in particular in the mixing matrices of gaugino couplings to fermions and sfermions. 

(hi) In the low energy effective theory, CP is an approximate symmetry. The Kobayashi- 
Maskawa phase is too small to account for ek which is explained, instead, by super- 
symmetric contributions. 
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(iv) The values of all other CP violating observables can be estimated and, in many cases, 
are drastically different from the Standard Model predictions. In particular, a^K s 
and a^vv are both much smaller than one. 

The mechanism that is responsible for the approximate nature of CP is basically the 
following. The full high energy theory has exact horizontal and CP symmetries. There are 
three relevant high energy scales: A#, where the horizontal symmetry is spontaneously 
broken; Acp, where CP is spontaneously broken; and Ap, where the information about 
these spontaneous breakings are communicated to the observable sector. There exists a 
hierarchy between these scales: Acp <C Ah <C Af (the hierarchy Acp <C Ah is naturally 
produced by the scalar potential) , so that in the low energy effective theory, the horizontal 
symmetry and the CP symmetry appear explicitly broken by small parameters: 

Ap 

This mechanism, while predicting phenomenology of CP violation that is very different 
from the Standard Model, also solves both the flavor and CP problems of supersymmetry. 

Our models use the Froggatt-Nielsen mechanism || to achieve the small breaking 
parameters. We employ supersymmetric abelian horizontal symmetries similarly to ||. 
The supersymmetric flavor problems are solved by the alignment mechanism |H|. As 
concerns CP violation in the supersymmetric framework, the idea of approximate CP has 
been discussed in refs. and spontaneous CP violation has been discussed in refs. 

||{l7],!|-[f|] . Our work is closely related to two of these works. In ref. J7|, models were 



constructed with spontaneous CP breaking and approximate CP in the low energy theory. 
However, while the mechanism of communicating the breaking in ref. is aimed to 
solve the strong CP problem and leads to a single low energy phase, our mechanism is 
aimed to solve the supersymmetric flavor problems and leads to a large number of low 
energy phases. Both the breaking mechanism and the communication mechanism are the 



same as in ref. [17]. The main new ingredient in our models is that, while the models 
of ref. JT7J have effectively CP breaking parameters of order one, the models presented 
in this work give small CP breaking and, therefore, a very different phenomenology of 
CP violation. Moreover, as the supersymmetric CP problem is solved partially by the 



approximate nature of CP, the required alignment is much less precise than in existing 
models. This situation gives more freedom in constructing the models and also allows for 
some different phenomenological signatures in FCNC processes. (We do not consider the 
strong CP problem in this work. Note, however, that the alignment models may solve this 
problem too [p~8[] . ) 

The structure of this paper is as follows. We first present two explicit models of 
approximate CP, one where the breaking parameter is intermediate, (9(0.04) (section 2) 
and the other where it is very small, (9(0.001) (section 3). The implications of these models 
for flavor changing neutral current processes are studied in section 4 and for CP violation 
in section 5. Section 6 clarifies an interesting point about holomorphic zeros, which bears 
consequences for rare K decays. Our conclusions are summarized in section 7. 



2. Model I 

Our first model employs a horizontal symmetry 

H = U(l) 1 xU(l) 2 . (2.1) 

The superfields of the supersymmetric standard model (SSM) carry the following H- 
charges: 

Qi(2,l) Q a (3,-1) Q 3 (0,0) di(4,-l) J 2 (-2,4) d 3 (l,l) (22) 
«i(5,-l) tZ 2 (-2,4) m 3 (0,0) d (-l,O) ^(0,0) { -> 

where Qi are the quark doublets, d{ and Ui are the down and up quark singlets, and fa are 
the Higgs doublet fields. In addition, we have three standard model singlet superfields: 

Si (-1,0), 32(0,-1), S 3 (-3,-l). (2.3) 

The horizontal symmetry is spontaneously broken when the three Si fields assume 
VEVs. The breaking scale is somewhat lower than a scale M where the information about 
this breaking is communicated to the SSM, presumably by heavy quarks in vector-like 
representations of the Standard Model 0. We will quantify all the small parameters as 
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powers of a small parameter A which we take to be of 0(0.2). Then, we take for the three 
VEVs 

Note that due to the U(l)i x £7(1)2 symmetry, we can always choose (Si) and (S2) to be 
real. However, (S3) is, in general, complex with a phase of 0(1). Then CP is spontaneously 
broken by £3. The hierarchy £3 <C ei, £2 and arg(£3) = 0(1) can be naturally induced, as 
explained below. 

The electroweak symmetry is spontaneously broken by the VEVs of (fid and (fi u , and 
we assume that 

The model is defined by the horizontal symmetry, the assigned horizontal charges and 
the hierarchy of VEVs. For most of our purposes, however, we need not consider the full 
high energy theory. It is sufficient to analyze the effective low energy theory, which is the 
SSM supplemented with selection rules that follow from the if-breaking pattern: 

(a) Terms in the superpotential that carry charge (m, n) under H with m, n > are 
suppressed by (9(A m+n ), while those with m < and/or n < are forbidden (due to 
the holomorphy of the superpotential 0). 

(b) Terms in the Kahler potential that carry charge (m, n) under H are suppressed by 
£>(AM+M). 

These selection rules allow us to estimate the various entries in the quark mass matrices 
M q and the squark mass-squared matrices M q2 . For each entry, we write the leading 
contribution and the subleading contribution if it is complex with respect to the leading 
one (namely, if it has a different £3-dependence). We do not write the coefficients of 0(1) 
which appear in each entry. For the quark mass matrices and for the off-diagonal blocks 
in the squark mass-squared matrices, we write the effective matrices after the rotations 
needed to bring the kinetic terms into their canonical form have been taken into account 
§. We get: 

ef 6ie\ t\e\ + £ 3 £i£ 2 

M d ~ ((fi d ) I e\e\ + e 3 e?e| ^4 ^ 4 + ^4 \ • (->•") 

e l e 2 + e 3 e 2 e l e 2 + e 3 e 2 e 2 



'el + e 3 (efel + efe 2 ) e| 
M u ~ (0.) I ef c| + e 3 e?ei eie§ e?e 2 + e 3 e 2 I , (2.7) 

t 5, i , ,2,2 ,2,4 I ,*, ,5 - ' 





/ 1 exel 

Mf 2 ~m 2 ei e| 1 efe 2 + e 3 e 2 ) , (2 A 

\efe 2 + e^ei e\t 2 + e\e\ 

( i 44 + 444 44+44\ 

M d A ~ rh 2 44 + e 3 el4 1 e\el + e 3 ef , (2.9) 

V 44 + es4 44 + 44 1 / 

/ 1 elel + egefel ef e 2 + e^e| \ 

~ ™ 2 e i e 2 + ^44 1 + e 3 eie| , (2.10) 

\efe 2 + e 3 e2 e 2 £ 2 £ 4 + £ * ei£ 5 1 J 

(Mf^- ~ m(M%, (2.11) 



We can also estimate the size of the bilinear \x and B terms: 

fi ~m(ei + 44 e 2), 

2 -2/ i * 2 \ 

m 12 ~m (ei + e 3 e 1 e 2 ). 



(2.12) 



Thus the horizontal symmetry solves the //-problem in the way suggested in ref. [|19 . 

From the mass matrices, we can further estimate the mixing angles in the CKM matrix 
and in the gaugino couplings to quarks and squarks. We denote the latter by where, 
for example, is the mixing matrix that describes the gluino couplings to left-handed 
down quarks and 'left-handed' down squarks. (The LR mixing angles are very small and 
we do not present them explicitly.) We write the estimates in terms of powers of A. For 
the CKM matrix, we find 

|K S |~A, |K & |~A 3 , |K b |~A 2 , (2.13) 
as required by direct measurements, and 

5 KM ~ A 2 . (2.14) 

For the gaugino couplings we find 

(K d L ) 12 ~ AV A4 , (i^) 13 ~AV A2 , (i^) 23 ~AV A4 , (2.15) 
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(Kl) 12 ~ A, (Kl) 13 ~ AV A , (K£) 23 ~ AV A , (2.16) 
(i^) 12 ~AV A2 , (Ki) 13 ~\*e a \ (K d R ) 23 ~ \*e iX \ (2.17) 
{K U R ) 12 ~ AV A \ (*T£) 13 ~ A 6 e lA2 , (K£) 23 ~ A 6 e^ 4 . (2.18) 



Note that in ( |2.15| )-( p.l8|) we omit coefficients of order one not only in the overall magnitude 



but also in the phases. 



Finally, we can estimate the relevant supersymmetric CP violating phases 20,21 : 



i B = axg(m? 2 //x) ~ A 6 , (2.19) 



= arg vgrn „ A 4 (2 20) 



while the corresponding <p A is negligible. 

Before concluding this section, we would like to show how a complex (S3) which is 
hierarchically smaller than (Si) and (S 2 ) can be achieved naturally. Let us add yet another 
Standard Model singlet field S^G, 2). The Si dependent terms in the superpotential are 

W ( S i) ~ ^^Sf + J^S 4 SfS 2 S 3 + cS 4 S 3 2 , (2.21) 

where a, 6, c are dimensionless numbers of 0(1). For (S4) = we have Fg 1 = i 7 ^ = = 
0, while Fs A = requires 



.6.2 , , , „,2 _ n _^ = ~ & ± ^ ~ 4aC 

e?e 2 2c 



ae^ + be\e 2 e 3 + ce\ = =^ = ^ . (2.22) 



We see that indeed |e 3 | ~ |ef €2 1 ~ A 4 and that for b 2 — 4ac < 0, e 3 is complex. (This 
mechanism for spontaneously breaking CP was first suggested in ref. @.) 



3. Model II 

Our second model employs a horizontal symmetry 



H = U(l) 1 xU(l) 2 xU(l) 3 . 
6 



(3.1) 



The SSM superfields carry the following if-charges: 

Qi(2,0,l) Q 2 (0,0,2) Q 3 (0,0,0) 4(0,5,0) 4(1,5,-2) 4(4,0,0) 
^(-2,6,0) u 2 (l, 1,0) %(0,0,0) d (-l,O,O) ^(0,0,0). 



(3.2) 



We have four standard model gauge singlet fields: 

Si (-1,0,0), S 2 (0,-1,0), S 3 (0,0,-1), 54(0,0,-4). (3.3) 
The orders of magnitude of the various S-VEVs are 

c ..- (Si) x c ... (S-l) x __ (S3) , (^4) , 4 fo A \ 

€l = ~m A ' e2 = IT ' 3 = ~W ' 4 = T ' ( } 

Due to the £7(1) 1 x Z7(l)2 x U(l)z symmetry, we can always choose (Si), (S2), (S3) real, 
but (S 4 ) is, in general, complex with a phase of 0(1). For the electroweak breaking VEVs, 
we take 

tan/9~l. (3.5) 

For the various quark and squark mass matrices, we get: 

/ eie|e 3 e 2 eie 3 (l + e§e|) ef e 3 (l + e§e 4 ) \ 
M d ~ <0 d > eie|e|(l + e 4 ) e| efe|(l + e 4 ) , (3.6) 

where M 3 d x ~ eie^e 2 + e§(l + e 4 + el)), 

4 e 3 efe 2 e 3 (l + e§e 4 ) efe 3 (l + e|e 4 ) 

M u ~ (0 U ) I ^£§(1 + 64) eie 2 e|(l + e 4 ) e§(l + e 4 ) | . (3.7) 

CiCn 6i€ 2 1 



/ 1 e?e 3 (l + €§€$) e 2 e3 (l + e 2 e4 )\ 

e 2 e 3 (l + e 2 e 4 ) 1 e 2 (l + e 4 ) , (3.8) 

Ve?e 3 (l + elea) e§(l + ej) 1 / 

eie§(l + ej) e\e\ \ 

Mf R ~ m 2 ( e ie |(l + 64) 1 ef e |e§(l + e 4 ) , (3.9) 

44 444(i + 4) 1 / 
1 e\4 44\ 

Mf R ~ m 2 I efe| 1 e ie2 , (3.10) 



e 1 e 2 eie 2 1 y 



(Ma) y ~m(M%. (3.11) 
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For the bilinear terms, we find 



\i ~mei(l + 63(64 + 64)), 

(3.12) 

m\ 2 ~m 2 ei(l + 63(64 + el)). 

For the CKM matrix, we find again magnitudes consistent with the measurements 
(namely, the same orders of magnitude as in ( |2.13| )) but the KM phase is smaller: 

8km ~ A 4 . (3.13) 

For the gaugino couplings we find 

(JtJ)ia ~ AV*", (i^) 13 ~AV A6 , (K d L ) 23 ~ AV A4 , (3.14) 

(K U L ) 12 ~ \e lX \ (KZ) 13 ~\ 3 e iX \ (K U L ) 23 ~ AV A \ (3.15) 

(K d R ) 12 ~ AV A \ (tf£) 13 ~ A 7 e lA4 , (i^) 23 ~ A 4 e lA \ (3.16) 

M) 12 ~AV A6 , (i^) 13 ~AV A< \ (K£) 23 ~A 2 . (3.17) 
The supersymmetric CP violating phases are: 

<Pb ~ A 8 , (3.18) 



tu ~ A 6 , (3.19) 



; .4 



while 0^ is negligible. 

The required hierarchy between the H and CP breaking scales is achieved by mini- 
mizing the Higgs potential for the four Si of eq. (|3.3|) and a fifth singlet field £5(0, 0,8). 
This would give 64 = 0(e 3 ) and complex. 



4. Flavor Changing Neutral Current Processes 

Generic supersymmetric models, with mass-squared differences between generations 
of 0(m 2 ) {fh is the supersymmetry breaking scale) and supersymmetric mixing angles of 
(9(1) give much too large contributions to various flavor changing neutral current (FCNC) 
processes such as Am^, Amo and Amg. There are various solutions to this problem: 

8 



a. All squark generations are equal at some high energy scale. This is the situation, for 
example, in models of gauge mediated supersymmetry breaking. 

b. The first two squark generations are degenerate due to a non-Abelian horizontal sym- 
metry. 

c. The first two squark generations are very heavy. 

d. Squarks are neither degenerate nor very heavy, but the mixing angles in the gaugino 
couplings to quarks and squarks are small. 

The last option arises naturally in models of Abelian horizontal symmetries of the 
type that we used in constructing our models. Indeed, one can easily see from eqs. fl2.15|) - 
( |2.18|) and ( |3.14|) -( [3.17| ) that there is no mixing angle of order one in our models; they are 
all suppressed by the selection rules of the horizontal symmetries. To understand whether 
the alignment in our models is precise enough to satisfy the phenomenological constraints 
and, in the case that it is, whether the supersymmetric contributions are significant in 
comparison to the Standard Model ones, we write down the constraints on the mixing 
angles (taken from ref. |22|J ) in terms of powers of A and then compare to the predictions 
of our two models. This is done in Table 1. (We define (Kij) = [(K^ijiK^jij] 1 / 2 .) 



Mixing Angle Process Bound Model I Model II 





ArriK 


A -A 2 


A 3 


A 3 


(K d R )i2 


ArriK 


A -A 2 


A 5 


A 3 


(Kf 2 ) 


ArriK 


A 3 


A 4 


A 3 


(K d L )is 


ArriB 


A 


A 3 


A 3 


(K d R )is 


ArriB 


A 


A 5 


A 7 


(Kf 3 ) 


Arris 


A 2 


A 4 


A 5 


{K u L )xi 


ArriD 


A 


A 


A 


TO) 12 


ArriD 


A 


A 4 


A 4 


m 2 ) 


Amp 


A 2 


A V2 


A 5 / 2 



Table 1. Supersymmetric mixing angles in our models and the phenomenological bounds on them. 
We learn the following points from the Table: 
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(i) The contributions to Amp that are proportional to [{K^)i2\ 2 saturate the experi- 
mental upper bound in both models. This is a generic feature of models of alignment 
related to the fact that in these models the Cabibbo mixing (|V^ S | ~ A) comes from the 
up sector. 

(ii) The contributions to Arris axe very small. In all alignment models, the stan- 
dard model amplitudes dominate ||. But while the supersymmetric contributions could 
be generically of (9(20%), the models constructed here provide an example where these 
contributions are below the percent level. 

(iii) The contributions to ArriK are of (9(10%) in model I and saturate the experi- 
mental value for model II. This is in contrast to all previous models of alignment where, 
to satisfy the Ek constraint, the supersymmetric contributions to ArriK were negligibly 
small. The large contribution comes in the two models from (Kf) 12(^^)12 • 

Before proceeding, we would like to make two comments: 

a. The contributions to other FCNC processes, such as ArriB a and b — > S7, are very 
small. The K + — > tv + uu decay is discussed separately below. 

b. The contributions from the (LR) blocks in the squark mass-squared matrices are much 
smaller than those coming from the mixing angles presented in Table 1. This is the 
reason why, even though we calculated them explicitly, we do not present them in 
Table 1. 

As concerns the rare K+ — > n + uu decay, the largest supersymmetric contribution in 
alignment models comes from (K^)i2, if it is as large as allowed by the ArriK constraint, 
(Kjj) 12 ~ A 2 [23]. In such a case, the supersymmetric contributions from penguin diagrams 
with chargino and u,c squarks are significant. In both our models, (1^)12 ~ A 3 , leading 
to supersymmetric contributions of 0(10%). While both the standard model and the 
supersymmetric amplitudes are real to a good approximation, so that there is maximal 
interference between the two, the relative sign is unknown so that the rate could be either 
enhanced or suppressed compared to the standard model. 

It is interesting that we are unable to construct a model where either (Kf)\2 or (-^^)i2 
are as large as allowed, namely C(A 2 ). This situation goes beyond the two specific models 
that we present here and seems generic to models with continuous Abelian horizontal 
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symmetries. The reason for that is explained in section 7. The situation is different in 
models of discrete Abelian symmetries. We actually constructed a model with a Z4 x Zg 
horizontal symmetry where (K^)^ ~ A 2 . The model is, however, quite complicated and 
its Higgs potential does not provide in a natural way the hierarchy between Ah and Acp 
(this seems a rather generic feature of models with discrete Abelian symmetries) , which is 
the reason that we do not present it here. 

The contribution to K + — > n + ui> from the (LR) sector is negligibly small. For the 
(LR) contributions to be significant, we need the off-diagonal terms in Mf iR to be of 
0(mtfh) ||24j| . while in our models they are much smaller than that, as can be seen from 
eqs. (pUD and (|3JTp. 



5. CP Violation 

Each of the two models that we have constructed has an approximate CP symmetry 
for the SSM. In model I, all CP violating phases are < 0(X 2 ) and in model II, they 
are < C(A 4 ). The resulting predictions for CP violating observables are then drastically 
different from the Standard Model. 

The first thing to note is that with 5km ~ A 2 or A 4 , it is impossible to account for 
Ek ~ 10 -3 by the Standard Model contributions. However, in both models, 

3[(i^) 12 (i^) 12 ]~A 10 , (5.1) 

which can account for ek from the supersymmetric gluino- mediated diagrams. 

The most dramatic consequences of the approximate CP symmetry concern the CP 
violating asymmetries that are expected to be large in the standard model. First, let us 
consider CP asymmetries in neutral B decays into final CP eigenstates. For the sake of 
definiteness, we consider a^K s - The supersymmetric contributions to the B — B mixing 
amplitude are, as mentioned above, negligible. Usually this leads to the conclusion that 
the standard model predictions for cl^k s remain valid. But this is definitely not the case 
in our framework. The fact that 5km is very small means that so will be cl^k s - Explicitly, 

<HKs_ = j 0{\ 2 ) Model I, , , 

^K S \°( X4 ) Modelll. 

11 



If a^Ks is measured to be in the Standard Model range, our models of approximate CP 
will be excluded. 

Concerning a nU i>, it was shown that a^^p = sin 2 9k, where 6k is the relative phase 
between the K—K mixing amplitude and the s — > dvv decay amplitude fl2~5" |. In our models, 



the standard model contributions dominate both the real part and the imaginary part of the 
decay amplitude. In model I, the standard model also dominates the real part of the mixing 
amplitude, while the supersymmetric contribution dominates the imaginary part. In model 
II, the two contributions to the real part are comparable, but supersymmetric diagrams 
dominate the imaginary part of the mixing amplitude. In either case, the approximate CP 
symmetry leads to a strong suppression of a^vv- 

Oiruu _ J C(A 4 ) Model I, 



„SM_ 



\ (9(A 8 j Model II. 



(5.3) 



We learn that if a^vv is measured in the foreseeable future, our models of approximate CP 
will be excluded. 

CP could play an interesting role in D — D mixing [26,27]. In measuring the time- 



dependent decay rate for D° — > K + tt~ and D° — > K~tt + , a term proportional to te~ rt that 
is different between the two CP-conjugate modes will appear if the relative phase between 
the D — D mixing amplitude and the c — > usd decay amplitude, 6d = arg X D o^ K + n - , 
is large. This is the case in previous alignment models. However, in our models, the 
relevant phase, that is arg[(iT^) 2 2 ], is very small (< (9(A 4 )), so the effect is probably 
unobservable. If D — D mixing is not observed within, say, one order of magnitude of the 
present experimental bound, then the existing alignment models are excluded. But if such 
mixing is observed and with large CP violation, then the alignment mechanism remains 
viable but not in combination with approximate CP. 

Finally, we discuss the electric dipole moment of the neutron d^- It was argued in ref. 
that in supersymmetric models without universality, namely when there is no super- 
GIM mechanism, there is a generic lower bound on the CP violating phases that contribute 
to djsf. This bound is of C(A 6 ) and leads to d^ > 10~ 28 e cm. This bound is about three 
orders of magnitude above the value in supersymmetric models with universality and may 
be within the reach of forthcoming experiments. Indeed, our models obey this bound and 
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predict a potentially observable d^. 

Our results concerning CP violation are summarized in Table 2. Note that, within 
the Standard Model, d^vv = C(sin 2 /3), which parametrically is of 0(1), but turns out 
to be numerically of 0(A) j2~g[j . d n is given in units of 10 -23 e cm, so that the present 
experimental bound is dw ?S A 2 . 



Process 


SM 


Model I 


Model II 




0(1) 


e>(A 2 ) 


e>(A 4 ) 




1 


C(A 4 ) 


C(A 8 ) 


9d 





< 0(X 4 ) 


C(A 4 ) 


djsr 





0(X 4 ) 


C(A 6 ) 



Table 2. CP violating observables in the SM and in our models. 

6. Lifting Holomorphic Zeros 

The Yukawa couplings, being part of the superpotential, are holomorphic in the H- 
breaking parameters. In particular, if all breaking parameters carry charges of the same 
sign under one of the horizontal £7(l)'s, then an entry in the Yukawa matrix that breaks 
H by a charge of the same sign vanishes 0. We call these vanishing Yukawa couplings 
'holomorphic zeros'. However, in the basis where these holomorphic zeros are exact, the 
kinetic terms are not canonical. When we normalize them back to a canonical form, the 
holomorphic zeros are lifted P,|29[|. 

Knowing the CKM mixing angles and the quark mass ratios allows us to guess a 'naive 
value' for each entry in the Yukawa matrices [Q. These are 

/A 4 A 3 A 3 \ /A 7 A 5 A 3 \ 

Y d ~A 3 tan/3 A 3 A 2 A 2 , r~ A 6 A 4 A 2 , (6.1) 
V A 1 1 J \X 4 A 2 1 / 

(where we have taken the high energy values m^/mt ~ A 3 and m c /mt ~ A 4 ). When one 
of these entries vanishes because of holomorphy, the rotation to bring the kinetic terms 
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back to the canonical normalization cannot lift this zero to its naive value ||. In fact, in 
our model, the entry in the effective Yukawa matrix (that is, in the basis where the kinetic 
terms are canonically normalized) is suppressed by at least A 2 compared to the naive value. 
Below we prove this statement. 

Before we provide the proof, we would like to point out a phenomenological con- 
sequence of this situation. The naive value ^12/^22 ~ sm$c ~ A is too large for the 
ArriK constraint (if there is no squark degeneracy). The way that alignment models solve 
this problem is by assuming that Yf 2 is a holomorphic zero. The Am^ constraint allows 
^12/^22 ~ If this bound were saturated, then the supersymmetric contribution to 
the K + — > n + uu decay could be significant ||23|| . What we learn, however, is that the 



maximal value that we can obtain for this ratio in our framework is Y^/l^l ~ (which 
is, indeed, realized in both models presented above). Therefore, in alignment models that 
employ continuous Abelian horizontal symmetries with small breaking parameters that are 
< O(X), the modification to the Standard Model prediction for K + — > n + uu is never large. 

To prove that a lifted holomorphic zero is suppressed by at least the square of the 
breaking parameters, let us consider a horizontal symmetry H — U(l) x x U(l) y broken by 
small parameters of order 

e x (-l,0)~A n -, e y (0,-l)~A^, (6.2) 
with two down-quark generations, 

Qi(a x ,a y ), di(b x ,b y ), Q 2 (cx,c y ), d 2 (d x ,d y ). (6.3) 
The kinetic terms for the (Qx, Q2) fields have coefficients of order 



I \n x \a x — c x \+n y \a y -c y \ 

\n x \a x —c x \+n y \a y —c v \ ^ 



(6.4) 



while those of the (di, d 2 ) fields will be 



]_ \n x \b x — d x \+n y \b v —dy\ 

\n x \b x — d x \+n y \b y — d y \ 



(6.5) 



We further assume, without loss of generality, that the Yukawa coupling Y^, which breaks 
H by charge (a x + d x , a y + d y ), is a holomorphic zero because a x + d x < 0, while the other 
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entries in Y d do not vanish: 

/ )Jl x (a x +b x )+n y (a v +b y ) q \ 

I X n x( c x+b x )+n y (cy+by) \n x (c x +d x )+n y (c y +d y ) I ■ 

(6.6) 

A straightforward calculation shows then that the effective Yf 2 is 
(Y l d 2 ) eS = A + B + C, 



_ \n x (a x +b x )+n y (a y +b y )+n x \b x — d x \+n y \b y — d y \ 
_ \n x (c x +d x )+n y (c y +d y )+n x \a x —c x \+n v \a y —c y \ 
(J _ yn x {b x +c x )+ny{by+Cy)+n x (\a x —c x \ + \b x —d x \)+n y {\a y —Cy\ + \b y — dy\) 



(6.7) 



This should be compared to the 'naive' value, which is 

(iia)nalve ~ \ n -( a '+ d ")+ n v( a v+ d v) . (6.8) 

Let us first examine A of eq. (|6.7|). We assume that d y > 6 y (otherwise the suppression is 
even stronger). From a x + b x > and a x + (i x < 0, we conclude that b x > d x . We learn 
that 

' rs, ^ 2 "™( b *- d ») < A 2 "". (6.9) 



naivi 

Second, we examine B of eq. (|6.7|). We assume that a y > c y . From c x + d x > 0, we 
conclude that c x > a x . We learn that 

TD 

/-y d \ — v ' ' 

\ 12 /naive 

As concerns C, it is easy to see that it is smaller than A and B. 

The final conclusion is then as follows. Suppose that a holomorphic zero is induced 
because the Yukawa coupling carries a negative charge under a symmetry U(l) x that is 
broken by a small parameter e x ~ X Ux . Then, the effective Yukawa coupling, that is the 
coupling in the basis where the kinetic terms are canonically normalized, obeys 



Pi 



<\ Zn *. (6.11) 



(Y q ) ■ ~ 

\ ijJ naive 

2 \ 



Since in all our models n x > 1, the suppression is at least by 0(A ). 
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7. Conclusions 

Supersymmetry allows for CP violating mechanisms that are dramatically different 
from the Standard Model. In particular, CP could be an approximate symmetry, with all 
CP violating phases very small, 

10" 3 < <Pcp < 1. 

In this work, we gave two examples of phenomenologically viable models where CP is 
broken by paramaters of order 0.04 or 0.001. 

The two models that we presented here are not unique. We use them to demonstrate 
how approximate CP can arise naturally and to explore the phenomenological signatures 
of approximate CP. The specific models should only be thought of as examples, but the 
underlying CP breaking mechanism and the resulting phenomenological implications are 
generic to this class of models. 

The fact that the Standard Model and the models of approximate CP are both viable 
at present is related to the fact that the mechanism of CP violation has not really been 
tested experimentally. The only measured CP violating observale, that is Ek, is small. Its 
smallness could be related to the 'accidental' smallness of CP violation for the first two 
quark generations, as is the case in the Standard Model, or to CP being an approximate 
symmetry, as is the case in the models discussed here. Future measurements, particularly 
of processes where the third generation plays a dominant role (such as a^K s or a„^), will 
easily distinguish between the two scenarios. While the Standard Model predicts large CP 
violating effects for these processes, approximate CP would suppress them too. 

The distinction between the Standard Model and Supersymmetry could also be made 
- though less easily - in measurements of CP violation in neutral D decays and of the 
electric dipole moments of the neutron. Here, the GIM mechanism of the Standard Model 
is so efficient that CP violating effects are unobservable in both cases. In contrast, the 
flavor breaking in supersymmetry might be much stronger, and then the approximate CP 
somewhat suppresses the effects but to a level which is perhaps still observable. 

Finally, we note that the predictions for FCNC processes are modified even for those 
processes where the supersymetric contribution is negligible. The reason is that the con- 

16 



straints on the CKM parameters are modified. Instead of the ek constraint, which is 
not relevant for the CKM parameters, we have rj ~ 0. The resulting modifications were 
recently analyzed in ref. |3(| . 
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